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Abstract 

We study the bottomonium spectrum in the nonrelativistic quark model with the coupled-channel 
effects. The mass shifts and valence bb component are evaluated to be rather large. We find that 
the hadronic loop effects can be partially absorbed into a reselection of the model parameters. No 
bottomonium state except T(5"'^S'o) and T(5^5i) with mass around 10890 MeV is found in the 
quark models both with and without coupled-channel effects, so we suggest that 1^(10890) is an 
exotic state beyond the quark model, if it is confirmed to be a new resonance. The predictions 
for the Xb(3-P) masses are consistent with the ATLAS measurements. If some new bottomonium- 
like states are observed at LHCb or SuperB in the future, we can determine whether they are 
conventional bottomonium or exotic states by comparing their masses with the mass spectrum 
predicted in our work. 

PACS numbers: 12.39.Jh, 12.40.Yx, 14.40.Pq, 14.40.Rt 
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I. INTRODUCTION 



In past years, the spectroscopy of heavy flavor quarkonium has seen great progress, partic- 
ularly the charmonium spectrum. Many charmonium-like states (such as X(3872), y(4260) 
and so on) with remarkable and unexpected properties have been reported. These exotic 
states present great challenges to our understanding of the structure of heavy flavor quarko- 
nium and quantum chromodynamics (QCD) at low energy, for a review, see Refs.[l-4|. 
On the other hand, many bottomonium states have been reported as well. In 2008, the 
spin-singlet pseudoscalar partner ?7f,(lS') was found by the Babar Collaboration with mass 
M = 9388.9li:3(stat) ± 2.7(syst) MeV The T{^Dj) was discovered in 2010 in the 
7r+7r-T(lS') final state with mass M = 10164.5 ± 0.8(stat) ± 0.5(syst) MeV |d]. In addition, 
the Babar Collaboration reported the P-wave spin-singlet hi,{lP) via its radiative decay into 
7?7fe(lS') with mass M = 9902 ± 4(stat) ± l(syst) MeV Q]- This state is confirmed by the 
Belle Collaboration j8[, and its mass is measured to be M = 9898.25 ± 1.06(stat)t|:[]^(syst) 
MeV. Meanwhile, the radial excitation state hh{2P) was also found by the Belle Collab- 
oration with mass M = 10259.76 ± 0.64(stat)t};o3(syst) MeV [Sj. Recently the ATLAS 
Collaboration has reported the discovery of the X6(3-P) state through reconstruction of the 
radiative decay modes of Xfe(3-P) — T(lS'j2S')7, and its mass barycenter is measured to be 
10.539 ± 0.004(stat.) ± 0.008(syst.) GeV ^. In particular, the Belle Collaboration has ob- 
served an enhancement in the production process e+e~ — )■ T(nS')7r"'"-7r~(n = 1, 2, 3) 10|. The 
fit using a Breit-Wigner resonance shape yields a peak mass of [10888. 4t2!g(stat) ±1.2(syst)] 
MeV and a width of [30.7j;7;o(stat) ± 3.1(syst)] MeV. In the following, we shall denote this 
state as 16(10890). Moreover, the Babar Collaboration measured the e+e~ — )■ bb cross section 



between 10.54 GeV and 11.20 GeV 



nj, the T(10860) and the T(11020) states, which are 



candidates of T(5S') and T(6S') respectively, were observed . Their masses and widths are fit- 
ted to be Mt(io860) = 10.876 ±0.002 GeV, rT(io860) = 43 ±4 MeV, Mt(iio2o) = 10.996 ±0.002 
GeV and rT(iio2o) = 37 ± 3 MeV, which are different from the previously measured val- 
ues. In particular, two charged narrow structures at 10610 MeV and 10650 MeV in the 
7r^T{nS){n = 1,2,3) and 7r^hi,{mP){m = 1,2) have been reported recently 12|. In sum- 
mary, the current experimental data indicate that there may be exotic bottomonium-like 
structures similar to the charmonium sector. Furthermore, LHCb has begun to run Il3l . 

n 

Belle will be updated to Belle II and a new SuperB factory will be built in Italy p^, we 
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expect that more heavy bottomonium states including the possible exotic extensions will be 
observed in the future. 

Motivated by the above exciting experimental progress in bb states, we shall carry out a 
careful, detailed study of bottomonium spectroscopy in this work, notably the poorly under- 
stood higher-mass bb levels. Thus, we can determine whether future observed bottomonium- 
like states could be accommodated as canonical bb states by comparing their masses with the 
mass spectrum predicted in this work. It is well-known that simple potential models, which 
incorporate a color coulomb term at short distances, a linear scalar confining term at large 
distances, and a Gaussian-smeared one-gluon exchange spin-spin hyperfine interactions, have 
been frequently used to describe both the charmonium and bottomonium spectrums. Gen- 
erally, the mixture between the quark model bb basis states and the two-meson continuum 
has been neglected in these models, which are called "quenched" quark models. 

The effects of the "unquenched quark model" including virtual hadronic loops have been 
studied extensively in the framework of the coupled-channel method 15-19|. The hadronic 
loop has turned out to be highly non-trivial, it can give rise to mass shifts to the bare 
hadron states and contribute continuum components to the physical hadron states. The 
possibility that loop effects may be responsible for the anomalously low masses of the new 
narrow charm-strange states D*q(2317) and 1^51(2460) has been suggested by several groups 
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-|23|. The hadronic loop in charmonium has been explored as well, and the mass shifts 
and continuum mixing due to loops oi D, D*, Ds and D* meson pairs have been studied 
extensively 24l-l3l|. Both the mass shifts and the two-meson continuum components of the 
physical charmonium states were found to be rather large. In particular, a J^^ = 1"^"*" state 
with mass about 3872 MeV could possibly be generated dynamically. 

Inspired by the large physical effects of hadronic loops in both the Dsj and charmo- 
nium states, we expect that the virtual hadronic loop should also play an important role 
in bottomonium spectroscopy. In this work, we shall study the bottomonium spectrum in 
detail, and take the hadronic loop effects into account. This paper is organized as follows. 
We present the framework of the coupled-channel analysis in section II. The non-relativistic 
potential model is outlined in section III. Section IV is devoted to the numerical results for 
the masses of the bottomonium states with and without hadronic loop effects, as well as 
phenomenological implications. We present our conclusions and discussion in section V. 
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II. FORMALISM OF COUPLED-CHANNEL ANALYSIS 




FIG. 1: Coupling of bb states to the BB mesons loop. 



In bottomonium, the process {bb) — {bn){nb) via light quark pair nn creation would 
induce the hadronic loop shown in FigUl where the initial bottomonium decays into inter- 
mediate virtual BI3 states and then reforms the original bottomonium state. Here B{B) 
denotes a general B{B) meson, it can be B{B), Bs{Bs), B*{B*) or B*{B*), the same conven- 
tion will be used henceforth without speci?cation. Since the open-flavor decay couplings of 
bottomonium states to two-body BB final states are large, the resulting loop effects should 
be important. This kind of virtual hadronic loop is universal, but it is not usually included 
in quark potential models and is only partially present in the quenched lattice QCD. The 
coupled-channel model is an appropriate framework for analyzing these hadronic loop effects 



15 



19| . In the simplest version of the coupled-channel model |25|, the full hadronic state is 
represented as 



with the normalization condition Y,a |cq 



^ + Ebb I d^P\XBB{pW = 1- denotes the 

bare confined bb states with the probability amplitude Cq, (pB (0_b) is the bn{bn) eigenstate 



describing the B{B) meson, and Xbb{p) is the wavefunction in the two-meson channel 
\4'b4'b)- The wavefunction obeys the equation 



n\'^) = M\^), n 




(2) 



where Hq is the Hamiltonian for the valence bb system, with the eigenstates determined 
by i^olV^a) = Ma\'4^a)- The Hamiltonian H^j^ acts between the constituents of B and B 
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separately, where the interactions between B and B are neglected. The continuum two- 
meson state \(f>B<pB) is the eigenstate of H^j^, 

Hbb\4>b4>b) = {Eb + EB)\(t)B(pB) - {rriB + ms + \'Pmi<Pm2) (3) 



BB 



where Eb = ym'^B + P^; -^s = + P^' ''^b and are the masses of B and B respec- 

tively, and fiBB = mB^m- i^ reduced mass of the two-meson system. Hi couples the 
bare state \ilja) with the two-body continuum Let us consider one bare state I'l/'o), 

the matrix element of Hi is of the following form: 

{<PBct>B\Hi\i^o)=hlB{v) (4) 

Substituting Eq.(l3]) and Eq.Q into Eq.([2]), we get the system of coupled equations for cg 
and Xi?B(p), 

CqMo + Y.j /i^b(p)xsb(p)c?'p = Mco 

BB-' 

{Eb + Eb)xbb{v) + Co/i^B(p) = ^Xbb(p) 
This coupled-channel equation can be solved straightforwardly, and we finally obtain the 
master equation 

M -M, + Y, IIbb{M) = Q. (6) 

BB 

Here H-BBiM) is the self-energy function for the hadronic loop induced by the intermediate 
states B and i?, it is explicitly given by 

IiBB{M) = I ^ '^^^^^1'! rf'p (7) 

^ J Eb + Eb - M -ie ^ 

Jsing the relation between the helicity amplitude /inn and the partial wave amplitude Aiis 



(5) 



32|, we have 

2 

B J^B 



2 











OO 2 771 771 

= r f dp- ^ V \Mls\' + ^^(^^^ E \Mls\') (8) 

J ^EB + EB-Mj^^ ' ^ M j^^ "^""^ J Eb+Ei,=m ^' 

Note that the imaginary part of the self-energy n5^(M) arises only if the initial hadron 
mass is above the intermediate BB threshold. Comparing with the two-body strong decay 



width shown in Eq. (lA.5p . it is obvious that the imaginary part is exactly equal to half of 
the decay width, if the decay is not forbidden kinematically. For the bottomonium state 
above the threshold, its mass is determined by the real part of the master equation Eq.([2]) 
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27, 



30|. The squared absolute value |cop is proportional to the probability that the 



physical energy eigenstate is in the bb configuration, and the bb component is given by 

In this work the effects of virtual hadronic loops will be considered in the above framework, 
and we shall sum over the contributions of intermediate loops from two stable S-wave mesons 
BB, BB*, B*B, B*B*, B^B,, B^B*, B*B, and B*B*, where antiparticles are indicated 
explicitly . 



III. NON-RELATIVISTIC POTENTIAL MODEL 

We use the standard non-relativistic potential model to describe the bare valence bb 
states. Its Hamiltonian is of the form 

Ho = ^ + Vcou{r) + Vcon{r) + C + Vsd{r) 
nib 

Vcouir) = Vconir) = ar, (10) 

O T 

where mi, is the bottom quark mass, Vcoui'^) is the well-known color-Coulomb force, and 
Konl'") denotes the linear confinement potential. To restore the hyperfine and fine structure 
of the bottom spectrum, one needs to introduce the Fermi-Breit relativistic corrections term 
VsdiT), which includes spin-spin, spin-orbit and tensor force. This is explicitly given by 

Vsd{r) = Vss{r) + VLs{r) + Vrir), (11) 

where Yssij"), Visij) and VTir) are spin-spin, spin-orbit and tensor operators respectively, 
and Vss is the contact hyperfine interaction, 

^''^'^ = ^-^^^^^-^'^^^ = ^ ■ 
where S;, is the spin of the bottom quark and is the spin of the anti-bottom quark. The 
Gaussian smearing of the hyperfine interaction is introduced here. 

Sir) = {^fe-^'^' (13) 



TT 
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The spin-orbit term is given by 

/ N r d^^ , N d^^ , L ■ S /4a, a.L-S ,^ 
V.s(r) = |3-l^„„(r) - = - -)^, (14) 

where S = S;, + is the total spin, and L is the relative angular momentum between h and 
h. Finally, the tensor term is 

where T^^ is the well-known tensor force operator, 

T,^ = 3(S,-r)(S5-r)-S5-S5 (16) 

The spin-dependent term Vsd is relativistically suppressed with respect to Vcou and Vcon, it 
is generally treated using leading order perturbation theory. 



IV. PREDICTIONS FOR THE BOTTOMONIUM SPECTRUM 



BB continuum, is an essentia 



well-established Pn model 
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The interaction Hamiltonian if/, which couples the bare valence hh with the two-body 

ement of this formalism. In this work, we shall use the 
] to describe the mixing between the bare hh state and the 
open bottom meson pair hq and g6. The ^Pq niodel assumes that the Okubo, Zweig, and 
lizuka (OZI) rule-allowed strong decay takes place via the creation of a quark- ant iquark pair 
with J^'^ = 0~^~^ from the vacuum. The qq pair production is described by the Hamiltonian, 

Hi = gY: /rf^x^^(x)^,(x) (17) 

where ipq is the Dirac quark field. Following the conventional calculating method in ^Pq 
model, one can then straightforwardly evaluate the valence-continuum coupling matrix el- 
ement /i^^ in Eq.(jlj). Here we will use simple harmonic oscillator (SHO) wavefunctions 
for the involved mesons, with a universal oscillator parameter /3. The SHO wavefunction 
enables analytical calculation of the transition amplitudes, and it turns out to be a good 
approximation. In Appendix A, we present the analytical ^Pq amplitudes M.ls for two 
S-wave final state channels. Since the oscillator parameters of the initial and final states 
are different in these analytical expressions, for an initial state of higher radial excited bot- 
tomonium, we can obtain the required transition amplitudes by simply taking derivatives, 



as is shown in Appendix B. Here we will take the typical values /3 = 0.4 GeV for numerical 



36-391, and it turned out 



calculations, this value was frequently adopted in the literature 
to be a reasonable zeroth-order approximation in coupled-channel calculations as well |27j . 
The parameter g has the form g = 2mq7, where is the constituent quark mass, and 
7 is the effective strength of pair creation. While for the creation of strange quarks, the 
effective strength 7s=(mg/ms)7 is used, following Ref. 25|, we shall take 7=0.322. The 
masses of constituent quarks are chosen to be m„ = = 0.33 GeV and = 0.55 GeV as 
usual. In the following, we first present the predictions for the bb spectrum in the nonrela- 
tivistic potential model of section III, then we include the coupled-channel effects, and the 
phenomenological implications are discussed. 



A. Bottomonium spectrum in conventional non-relativistic potential model 

As usual, we can determine the bare bb mass spectrum by solving the following 
Schrodinger equation numerically, 

- —^Uni{r) + [V{r) + ^l±^]uniir) = E^mUr) (18) 

where E^i is the energy eigenvalue, the potential V{r) is the leading order one with V{r) = 
Vcou{r) + Vcon{r) + C . The wavefunction '^nim of the system is closely related to Uniir) by 

^nim(r, e, <P) = R^i{r)Yi^{e, cp) and u„;(r) = r/2„i(r) (19) 

The mass splitting within the multiplets is determined by the spin-dependent term Vsd{r), 
which is taken to be a perturbation. As a result, we need to calculate the expectation value 
of Vsd{r) between the leading order wavefunction "^nim, in order to get the bottomonium 
mass fine splitting. The only term which gives a nonvanishing contribution to n ^Sq — n ^Si 
fine splitting is the spin-spin term. Since the expectation value of S;, ■ in the ^5*0 and ^5*1 
states is —3/4 and 1/4 respectively, the bare masses of S-wave states are 

M^{v}So) = E^o - l^^CSir)) (20) 

Main's,) = E^o + l^CS{r)) (21) 



where 



{6{r)) = {—r e-'^\l,{r)dr (22) 

^y^T Jo 



For the n ^Pi — n ^Pj, n — n ^Dj, n ^F^ — n ^Fj cases, the spin-orbital and tensor terms 
contribute as welL The corresponding bare masses are given by 



Mo(n^Pi) 


= Enl 


4 ^^^'^^ 






(23) 


Moin^Pj) 


= Enl 


+ 4 ,,nl + 


— [Aja,{—) H 


-Ejo{-)) 


(24) 


Mo{n^D2) 


= En2 


6 6ZTxas 

4 9mg ^^^'^^ 






(25) 




= En2 


"1 oo 

+ 4 + 


1 1 

— ) - 


-1 


(26) 


Moin'F^) 


= En3 


6 dZnas .7., 
4 9mg ^^^'^^ 






(27) 


Moin'Fj) 


= En3 


1 327ras ~ , , , 
+ 4 9,n| <^''')> + 


1 / ,1 


-Fja{-)) 


(28) 




= Eni 


4 9mg ^^^'^^ 






(29) 


Moin'Gj) 


= En4 


1 32-7ras ~ , , . 
+ 4 9™g + 




(30) 



Here (...) denotes the expectation value, which is defined in the same way as that in Eq.f l22|) . 
The coefficients Aj, Bj, Cj, Dj, Ej, Fj, Gj and Hj for each case are listed in Table [B 
Note that the tensor force Tf,i could lead to the so-called S— D mixing. However, the mass 
splitting induced by the S— D mixing is a second-order perturbation effect of the hyperfine 
interactions Vsd, thus its contribution is rather small and hence is neglected here. We find 
that the numerical value {S{r)) is extremely small except for the S*— wave states, since it 
would be proportional to the zero point value of the wavefunction, i.e., it is exactly zero if 
we don't smear the hj^erfine interaction. As a result, the following sum rules are satisfied 
quite well, 

Moin'Pi) - {5Mo{n^P2) + 3Mo{n^Pi) + Mo{n^Po))/9 (31) 
Mo{n'D2) ~ {7Mo{n^D3) + 5Mo{n^D2) + 3Mo{n''Di))/15 (32) 
Moin'Fs) ~ (9Mo(n3F4) + rMoin^F^) + 5Mo{n^F2))/21 (33) 
Moin^G^) - (UMoin^G^) + 9Mo{n^G4) + 7Mo{n^G3))/27 (34) 

For experimental input we use the masses of the 14 reasonably well-established bb states, 
which are given in Table UTTl Here we don't include the T{5S) and T{6S) candidates T(10860) 
and T (11020), since their masses measured by the Babar Collaboration are different from 
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J 


Aj 


Bj 


Cj 


Dj 


Ej 


Fj 


Gj 


Hj 





-16/3 


1 














1 


-4/3 


1/2 


-20/3 


3/2 










2 


28/15 


-1/2 


-4/3 


1/2 


-128/15 


2 






3 






80/21 


-1 


-4/3 


1/2 


-220/21 


5/2 


4 










52/9 


-3/2 


-4/3 


1/2 


5 














256/33 


-2 



TABLE I: Numerical values for the A j, Bj,Cj, D j, Ej, Fj, G j and Hj coefficients. 



the Particle Data Group (PDG) averages 40|. The parameters that follow from fitting these 
masses are a,=0.3840, a=0.9155GeV/fm, C = ^.7825 GeV, mfe=5.19 GeV and k=2.3 



41| , the effective strong coupling constant 



GeV. In the classical Godfrey-Isgur quark model 
as at the bottomonium scale and the string tension a are determined to be about 0.25 
and 0.18 GeV^ respectively. Clearly the fitting value for a is approximately the same as 



that in 4l|, while is found to be somewhat larger than that of 4l|. Note that is 
dissociated from the effective strong coupling constant determined by the hadronic width 
of the quarkonium in conventional potential models, it is only a purely phenomenological 
strength parameter of the short-range potential 19|. Given these values, we can predict the 
masses of the currently unknown bb states. The predicted spectrum is shown in Table IIIIl 
and we see that the mass sum rules in Eq. (l3Tl) -Eq. flMl) are really satisfied rather well. 



B. Bottomonium spectrum with coupled-channel effects 

Following the formalism presented in section II, we shall take the coupled-channel effects 
into account. By performing a fit to the 14 established experimental states given in Table 
TTTl the best values of the parameters are determined to be 0^=0.418, a=0.818GeV/fm, 
C = —0.62376 GeV, mh=5. 18 GeV and k=2.85 GeV. Note that the value of as here is 
larger than that in Ref. 411 , while a is smaller than the fitting value of 411 • We can now 
straightforwardly evaluate the bare state masses, the mass shifts due to BB loops and the 
predictions for the bottomonium masses with coupled-channel effects included. The results 
are given in Table IIVI and Table |Vl where the bb component is presented as well. It can be 
seen that the mass sum rules in Eq. ( !3T|) -Eq. ( l34l) remain approximately intact. In order 
to clearly see the predictions for the bottomonium spectrum, we further plot the predicted 
masses in Fig. [21 
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FIG. 2: Predicted and observed spectrums of bottomonium states. The solid lines are experi- 
mentally measured data, and the dashed lines are predicted masses with coupled-channel effects 
included. Various BB thresholds are also shown. 

From Table IIVI and Table El we can see that the mass shifts are predicted to be close 
to each other, they are of order 100 MeV, although individual BB loops make different 
contributions to the mass shift of each state. Moreover, we see that the total BB components 
are rather large, they mostly scatter in the range of 0.1 ~ 0.25. We conclude that the results 
for the hadronic loop contributions to the bottomonium states are consistent with the general 
loop theorem derived in Ref. 27l|. We note that the mass shifts in the charmonium sector 
are predicted to be around 200 MeV (25| (a larger value of about 500 MeV was suggested 



m 



27|), and the two-meson continuum components can be as large as 0.5. Therefore the 



hadronic loop effects in bottomonium are much smaller than those in the charmonium sector. 

In order to compare the mass spectrum in the non-relativistic potential model, we show 
the predicted mass with coupled-channel effects in Table IIIII as well. For the 16 observed 
states, both predictions agree well with the measured masses except T(4^5'i), T(5'^5'i) and 
T(6'^S'i). We see that the masses of T(4^S'i) and T(6^S'i) are predicted to be rather close 
to the observation after the coupled-channel effects are included, but the mass of T(5 ^Si) is 
smaller than its measured value about 51 MeV. Whereas the non-relativistic potential model 
without hadronic loop predicts the T(5'^5'i) mass successfully, the departures of predictions 
from observation are 41.6 MeV and 57.4 MeV respectively for the T(4^S'i) and T(6'^5'i) 
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states. The current experimental data indicates that coupled- channel effects can improve 
the agreement with observations, although it is not as satisfactory for the T(5 ^Si) state. But 
this is not the whole story, because the recently measured masses of T(5 ^Si) and T(6 ^5*1) by 
the Babar Collaboration are different from the world averages. More precise measurements 
of the T(5 '^S'l) and T(6 ^Si) states would be helpful in understanding the unquenched effects 
induced by hadronic loops in the bottomonium spectrum. Regarding the mass barycenter 
of the Xfe(3-P) states, it is predicted to be 10523.3 MeV and 10537.9 MeV respectively for 
the included and not included coupled-channel effects. Both predictions are consistent with 
the ATLAS measurement, although the former is slightly smaller. 

Inspecting the mass predictions shown in Table UTTl we can't find a state with mass around 
10890 MeV except T(5^S'o) and T(5'^5'i), even if coupled-channel effects are considered. 
Therefore if 10890) is confirmed to be a new resonance by future experiments, it should 
not be a canonical bottomonium state. Many interpretations have been proposed so far. 



A possible explanation is the existence of a tetraquark state [42|, |43|. Another 

explanation is a bb counterpart to the y(4260) state, which may overlap with T(5S') 44 1. 
Other interpretations such as final state interactions and so forth have been suggested as 



well 



45j. In the same manner, if some new bottomonium-like state is observed by LHCb 



or SuperB in the future, we can determine whether the state could be accommodated as 
a conventional quark model state by comparing its mass with our predictions, and we can 
determine its assignment if it is. 

Mixing between two bare bb states could arise through the hadron loop, provided that 
both bottomonium states could couple to the same intermediate BB state. Similar to the 
S— D mixing, this kind of mixing would introduce corrections to the mass of the physical 
bottomonium. However, its contributions are of higher order in the valence-continuum cou- 
pling Hamiltonian Hj with respect to the loop contributions discussed above. Moreover, as 
is stated by the loop theorem of Ref . , if the mass difference between various intermedi- 
ate loop mesons is neglected, the mixing amplitude between two valence bb states vanishes 
unless both the orbital angular momentum and the spin of the two states are the same. 
Concrete numerical calculations show that this mixing effect is really quite small [^^j, and 
the same results are found in Ref. 25|]. As a result, we have not considered this effect in 



the present work. 

Finally, we note that although the bottomonium mass spectrums predicted in models with 
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and without coupled-channel effects are not drastically different from each other, i.e., part of 
the hadronic loop effects can be absorbed into the redefinition of the model parameters, the 
underlying physics is different. In the scenario with coupled-channel effects, there are sizable 
BB components in the physical bottomonium states. As a result, the production and decay 
of bottomonium would be different from those predicted in the quenched quark model. The 
coupled-channel effects in hadronic transitions of bottomonium have been studied in Ref. 



461], and it was found that the inclusion of the coupled-channel effects can really improve 



the theory of hadronic transitions. The mixture of BB continuum may also be important 

)bservations, since the same has turned out to be true 
|. This topic deserves much laborious and complex 



in understanding some anomalous o 
in the charmonium sector y, Q, 1471, 



work and is beyond the scope of the present work. 



V. CONCLUSIONS AND DISCUSSIONS 



Motivated by the recent experimental progress on the bottomonium spectrum,, we spec- 
ulate that some bottomonium-like states may be observed in the future similar to the char- 
monium sector, particularly with the running of LHCb and SuperB. Because the coupled- 
channel effects are very important in understanding the nature of the newly observed 
charmonium-like states, such as X(3872), the same is expected to be true in the bottomo- 
nium sector. In this work, we investigate the hadronic loops effects in the bottomonium 
spectrum. The coupling between the valence bb states and the two-meson BB continuum 
is described in terms of ^Pq model. The mass shifts and the bb component of the physical 
bottomonium are calculated in detail. We find that the mass shifts for all states are similar 
to each other, they are around 100 MeV, although the contributions of individual loops are 
different for each state. The two-meson continuum BB components are found to be rather 
large as well. The hadronic loop effects in bottomonium turn out to be smaller than the ones 
in the charmonium sector. Moreover, we evaluate the mass spectrum in the conventional 
constituent quark model, where coupled-channel effects are not taken into account. We find 
that the hadronic loop effects can be partially absorbed into a reselection of the model pa- 
rameters. Since the potential models both with and without coupled-channel effects don't 
predict a state with mass around 10890 MeV except for T(5 ^Sq) and T(5 ^^i), we conclude 
that 1^(10890) should be an exotic state beyond the quark model, if it is confirmed to be a 
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new resonance. Our prediction for the mass of Xbi^P) states is consistent with the recent 
ATLAS measurement. Moreover, with the mass spectrum predicted in our work, we can de- 
termine whether a new bottomonium-hke state observed in the future can be accommodated 
as a canonical bottomonium, and we can determine its assignment if it is. 
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Appendix A: ^Pq transition amplitudes 



The derivation of the Pq matrix elements has been discussed in detail in Ref. [36 
Starting from the pair production Hamiltonian Hj given in Eq. flTTl) . one can straightfor- 
wardly evaluate the Hj matrix element hfi for the transition A{bb) — ?■ B[qb) + C{bq) in 
terms of overlap integrals in flavor, spin and spatial spaces, 

hfi ={BC\Hi\A)a + {BC\Hj\A), 

I signatureifl^^ flavorij^I spin+spaceis^ ~l~ signatureip^I flavorip^I spin+spacei]-)) (-^-l) 

Here a and h represent two decay diagrams in which the produced quark goes into meson B 
and meson C respectively. In the present work, only the first diagram is allowed. Therefore 
we have the flavor factor Ifiavor{o,) = 1 and Ifiavor{b) = 0, as is listed in Table [ITl The 
spin-space part of hfi is explicitly given by 

Ispin+space{a) = J d^k ^u^LaM^^ (k - Vb)^<1^l^^,j^^ (k - rPe) 



where r = and denotes created quark mass m„, or m^. '^uaLaMl^ is the 

wavefunction of the initial meson A in momentum space, and ubLbMl^ ^"^^ ^ucLcMl^ 
the wavefunctions of the final state mesons B and C respectively. Taking into account the 
phase space, we get the differential decay rate 

dn Ma ' ^ ' ^ ^ 

where P is the momentum of the final state mesons in the rest frame of meson A 



P = ^[Ml - {Mb + McY][Ml - {Mb - McY]/{2Ma). (A.4) 
To compare with the experiments, we transform the amplitude hfi into the partial wave 



amplitude ls by the recoupling calculation 



32j, then the decay width is 



T{A^B + C)= 27r— ^ ^ \M,s?. (A.5) 



Since we neglect mass splitting within the same isospin multiplet, to sum over all channels, 
one should multiply the mass shift due to a specific hadronic loop by the flavor factor J-" 
which is listed in Table HIl 
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Generic Decay 


Example 


I flavor («) 


Iflavor{b) 


T 


Y 

Jv - 


_^ R R 
— ? JDD 


Y 


— r D ^ JD 


1 





2 


X - 


^ B*B 


X - 


S*+ + B~ 


1 





4 


X - 


> B*B* 


X - 


S*+ + B*- 


1 





2 


X - 


> BsBs 


X 


^ + m 


1 





1 


X - 


> BlBs 


X 


Bf + B^ 


1 





2 


X - 


y D* R* 


X - 


Bf + Bf 


1 





1 



TABLE II: Relevant flavor weight factors for bottomonium decay, where \X) = \bb). 



We take all spatial wavefunctions to be simple harmonic oscillator forms with Pa being 
the oscillator parameter of the inital meson A and (5b = Pc for the final state mesons B and 
C. It turn out that the transition amplitude J^ls is proportional to an overall Gaussian 
factor, it can be expressed as 

Mls = ^e"^#^A5 (A.6) 



where r = — — , m„ denotes created quark mass rriu, rrid or m^, and P is the momentum 



of the final state mesons in the rest frame of meson A. The analytical expressions of the 
amplitudes for the decays into two stable S*— wave final states are listed in the following. 



37|. The mass difference 



We note that our expressions are different from the results in Ref. 
between the created quark and bottom quark and two oscillator parameters corresponding 
to initial and final states are considered in our expressions, and our results are more general 



than those in Ref. [37|. In the limit oi Pa = Pb = P and r = 1/2, the amplitudes presented 



below coincide with those of Ref. 
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Appendix B: Recursion relations between wavefunctions 



In Appendix A, we present the amplitudes for the ground bottomonium states decaying 
into two S— wave final states. If the initial state is the radial excited bottomonium, the 
corresponding amplitudes can certainly be derived in the same way. However, it is inter- 
esting to notice that the radially excited wavefunctions can be related to the lowest radial 
wavefunctions by differentiation, 
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/2' 1 2r? _ 2 8 

^^^^^P) = 5! ^^^-^ + ^)>^iA/.(p)e-^ = y 5/3^*ip(P) (A.67) 



v[/3p(p) = - ^ + 35)v[/ip(p) = 4^(5 + 2/^4 + 2/32^)v[/ip(p) (A.68) 



2^70' /34 /32 ^ ' ' 8(5 ' 8(5 



/8' 1 12 8 

*2d(p) = ^fj,-^Pi-^ + ^)>'2MjP)e"^ = ]jf-^1'Mp) (A.69) 



1 2n^ \/2 (9 

*2f(p) = ^(-9 + ^)*if(p) = ^/5^*iHp) (A.71) 

The differential operators depend only on /3, hence they can be pulled out of the integration 
in Eg. flA.2p . Therefore the amplitudes for radially excited meson decays can be found 
by applying the above differential operators to the amplitudes listed in Appendix A. For 
example, 

MLsi2 ^Si ^1 So So) = ^Pa^MlsC^^ So +' So) (A.72) 

VO op A 

We have checked that the amplitudes obtained with this method are exactly the same as 
those that result from performing the overlap integral straightforwardly. 
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TABLE III: Th e sp ectrum of the bottomonium states, where M^x is the PDG average for the 



measured mass 40(|. M„p denotes the prediction for the mass of the bb state in the conventional 
non-relativistic potential model presented in section III. Alth is the theoretical prediction after the 
coupled-channel effects are taken into account. All the masses are in units of megaelectronvolts 
(MeV). 
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0.908 





2.528 





1.727 


1.017 








2.332 


BtBs 


5.293 


3.519 


7.102 





6.011 


4.952 


3.879 


6.763 


6.346 





BlBl 


5.281 


5.971 


6.797 


10.795 


7.738 


7.456 


6.488 


7.941 


6.003 


9.510 


oM 


55.498 


58.246 


85.741 


81.826 


84.778 


87.333 


68.037 


97.692 


78.781 


75.047 


Mr, 


9447.3 


9518.6 


10001.2 


9957.1 


9991.6 


10016.9 


1 nriQzL 9 


1 n9zL9 


luoo 1 .y 




Mth 


9391.8 


9460.3 


9915.5 


9875.3 


9906.8 


9929.6 


10026.2 


10144.6 


1 no m 1 

10259.1 


10227.9 




9390.9 


9460.3 


9899.87 9859.44 


9892.78 


9912.21 


10023.26 


10164.5 


10259.8 


10232.5 


Phi 


0.877 


0.959 


0.910 


0.917 


0.912 


0.907 


0.923 


0.870 


0.885 


0.894 




X6i(2P) Xb2i2P) 


T{3S) 


T(45) 


T(5S') 




Vb{2S) 




%(4S) 


%(55') 


BB 





8.165 


5.441 


6.548 


8 Q63 


9 884 














B*B 


31.697 


22.568 


19.828 


24.377 


28.159 


21.871 


28.813 


28.9541 


33.6972 


39.649 


B*B* 


33.951 


37.106 


31.938 


34.353 


35.684 


48.760 


26.783 


26.720 


28.662 


29.955 


BgBs 





1.554 


0.992 


1.014 


1.048 


1.145 














BtBs 


5.410 


4.400 


3.756 


3.770 


3.952 


4.064 


5.694 


5.533 


5.573 


5.809 




6.804 


6.611 


6.215 


6.178 


6.424 


6.729 


5.447 


5.238 


5.224 


5.410 


6M 


77.863 


80.403 


68.174 


76.296 


84.229 


85.455 


66.215 


66.446 


73.155 


80.822 


IVlQ 


10330.3 


10350.5 


10420.1 


10679.0 


10904.2 


11108.1 


10071.26 


10404.4 


10666.4 


10893.4 


Mth 


10252.4 


10270.1 


10351.9 


10602.7 


10819.9 


11022.6 


10004.9 


10337.9 


10593.2 


10812.6 




10255.46 10268.65 10355.2 


10579.4 


10865.0 


11019.0 












0.887 


0.881 


0.893 


0.638 


0.791 


0.891 


0.926 


0.898 


0.776 


0.780 




77b(65) 


hiSP) 


X6o(3P) Xw(3P) Xb2{3P) T{1^D2) T{1^D,) T{l^Ds) T{2^D2) T{2^D^) 


BB 








17.469 





9.065 





13.077 


10.995 





13.187 


B*B 


36.924 


36.078 





33.683 


23.284 


43.057 


11.968 


27.296 


40.150 


11.383 


B*B* 


41.997 


31.981 


46.976 


33.449 


37.752 


40.082 


57.205 


45.513 


36.199 


51.163 


Bg Bs 








2.291 





1.499 





1.817 


2.088 





1.573 


BtBs 


5.980 


6.077 





5.207 


4.195 


7.575 


1.703 


5.272 


6.657 


1.445 




5.611 


5.688 


8.980 


6.423 


6.282 


7.175 


11.023 


7.531 


6.234 


9.663 


6M 


90.511 


79.824 


75.716 


78.762 


82.077 


97.890 


96.793 


98.696 


89.239 


88.412 


Mo 


11098.5 


10603.0 


10571.6 


10596.1 


10614.5 


10243.4 


10234.9 


10248.0 


10517.1 


10508.8 


Mth 


11008.0 


10523.2 


10495.9 


10517.3 


10532.4 


10145.5 


10138.1 


10149.3 


10427.9 


10420.4 


Mex 

Pbb 


0.883 


0.831 


0.843 


0.835 


0.816 


0.870 


0.872 


0.869 


0.836 


0.835 



TABLE IV: Mass shifts (in MeV) and bh component P^^i of bottomonium states due to individual 
BB loops. Mex is the PDG average for the measured mass, Mo denotes the bare mass, and Mth 
denotes the predicted masses with coupled-channel effects considered. 5M denotes the total mass 
shift. For simphcity, we have abbreviated the BB hadronic loop as ''BB" , BB* + B*B as ''BB*''\ 
and so forth. 
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T{3^D3) Til^Fs) T(i3F2) TCl^Fg) 


BB 





1 n 460 





1 3 300 





1 1 995 





1 3 733 





B*B 


38.298 


25.121 


43.392 


13.795 


43.215 


24.146 


46.402 


16.135 


44.945 


B*B* 


37.940 


41.427 


41.000 


54.099 


41.229 


45.943 


42.246 


58.809 


43.726 


BgBs 





1.879 





1.564 





1.841 





1.631 





B*sBs 


5.892 


4.690 


6.396 


1.397 


5.664 


4.531 


7.435 


2.001 


6.831 


B,B, 


f\ OA C 

o.y4o 


6.465 


5.921 


9.211 


6.601 


6.181 


6.965 


10.636 


7.542 


OlVl 


89.078 


90.041 


96.710 


93.367 


96.709 


93.936 


103.048 


102.945 


103.043 


Mr) 


10515.9 


10521.5 


10752.4 


10744.3 


10751.2 


10756.7 


10425.1 


10422.5 


10425.0 


Mil, 
i\ J 


10426.8 


10431.4 


10655.7 


10650.9 


10654.5 


10662.8 


10322.1 


10319.5 


10322.0 




0.836 


0.835 


0.922 


0.645 


0.763 


0.588 


0.829 


0.827 


0.829 




Til^Fi) T{2^F3) T{2^F2) T{2^F^) T{2^Fi) T{l^Gi) Tfl^Gg) T{I^Ga) Tfl^Gg) 


BB 


12.378 





18.347 





13.133 





15.860 





15.784 


B*B 


28.246 


47.250 


17.108 


46.383 


27.719 


49.311 


19.559 


48.645 


31.174 


B*B* 


48.263 


40.546 


55.203 


41.679 


47.068 


43.356 


58.613 


44.229 


50.179 


BgBs 


2.194 





1.423 





2.066 





1.520 





1.517 


BtBs 


5.142 


6.757 


1.699 


6.140 


4.789 


7.005 


2.050 


6.565 


3.272 


B* B* 


7.175 


6.257 


9.705 


6.844 


6.288 


6.475 


9.815 


6.900 


6.608 


6M 


103.398 


100.810 


103.484 


101.047 


101.064 


106.147 


107.417 


106.338 


108.533 


Mo 


10426.5 


10667.2 


10664.4 


10667.1 


10668.9 


10579.4 


10579.2 


10579.7 


10579.3 


Mth 


10323.1 


10566.4 


10560.9 


10566.1 


10567.9 


10473.3 


10471.8 


10473.4 


10470.8 


Mex 


0.830 


0.751 


0.637 


0.748 


0.749 


0.775 


0.762 


0.773 


0.733 



TABLE V: The continuing of Table HVl 



29 



